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OKlGiWAL PAGE 13 

OF POOR QUALITY 


OPTIMAL SYMMETRIC FLIGHT WITH AN INTERMEDIATE 

VEHICLE MODEL 

by 


P. K. A. Menon 
(ABSTRACT) 

Optimal flight in the vertical plane with a vehicle model 
intermediate in complexity between the point-mass and energy 
models is studied. Flight-path angle takes on the r6le of a 
control variable. Range-open problems fearure subarcs of 
vertical flight and singular subarcs as previously studied. 

The class of alti tude-speed-range-time optimization 
problems with fuel expenditure unspecified is investigated 
and some interesting phenomena uncovered. The maximum-lift- 
to-drag glide appears as part of the family, final-time- 
open, with appropriate initial and rerminal transient 
maneuvers. A family of climb-range paths appears for thrust 
exceeding level-flight drag, some members exhibiting 
oscillations. Oscillatory paths generally fail the Jacobi 
test for durations exceeding a period and furnish a minimum 
only for shorr-duration problems. 

Minimizing paths of long duration follow a cerrain 
corridor in the V-h chart. The features of the family 



sharpen for the special case of thrust and drag independent 
of altitude, and considerable analytical attention is 
accorded to this for the insight it provides to the more 
general model. 

The problem of "steepest climb" is found to be ill- 
posed with the vehicle model under consideration, straight- 
vertically-upward maneuver sequences being furnished by a 
family of paths alternating between upward and downward 
vertical flight and including a limiting "chattering" 


member . 
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Chapter I 
INTRODUCTION 

There has been interest from the begining of optimal flight 
studies in approximations featuring simplified vehicle 
models (Refs 1 through 8). The motivation in these studies 
stems from the considerable difficulties encountered in 
optimizing atmospheric flight trajectories with a detailed 
system model. In general, straightforward application of the 
maximum principle to an "exact” model of an aircraft leads 
to a nonlinear two-point-boundary-value problem involving 
several unknown parameters, with attendant difficulties in 
obtaining solutions. Further, a combination of three control 
variables, vis., lift coefficient, bank angle, and throttle 
and the imposition of various constraints considerably adds 
to the complexity of the problem. Nevertheless, it is 
currently feasible to compute control histories, i.e., open- 
loop controls, for optimal control of sophisticated 
mathematical models of flight vehicles. The favored approach 
involves reformulating the problem as one of nonlinear 
programming by parameterizing the control history (Ref 9). 
Such direct methods are computationally less sensitive than 
the indirect methods of solution via the state-Euler system 
with split boundary conditions. With the availability of 
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improved numerical techniques such as multiple shooting 
(Refs 10, 11 and 12), indirect methods are again begining to 
appear attractive for solving aircraft trajectory 
optimization problems (Ref 13). 

While numerical optimization techniques are of 
considerable value, they are expensive to use and with the 
current state-of-the-art in computing technology (Ref 14), 
virtually impossible to implement on-board an aircraft. 
Moreover, the results obtained are restricted to a 
particular set of boundary conditions and when these are 
altered, the numerical exercise must be carried out all over 
again. It may be added that these techniques are ill-suited 
for use in preliminary design stages where it is desirable 
to have the capability to assess the effect of design 
changes rapidly and economically. 

Experience from flight tests as well as comparisons 
between various solutions of the optimal-control problem 
often reveals that the improvement in performance is minimal 
when the exact optimal trajectory is compared with a 
suboptimal one obtained from simpler analysis (Ref 15). 
Specifically, for the purposes of flight-path optimization 
for transport aircraft that do not generally indulge in 
violent maneuvers, reduced-order models have been found 
adequate (Ref 16). The simplified analysis has the added 
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advantage that the resulting solution is close to optimal 
and can sometimes be corrected for the effects of neglected 
dynamics. The use of simplified approximations is not free 
of difficulties, however. For instance, in Ref 17., optimal 
control resulting from the introduction of various 
approximations into the equations of motion was compared 
with the optimal control based on an "exact" set of motion 
equations and it was found that, in some cases, the 
resulting controls violated the assumptions made. It is 
prudent, therefore, to exercise caution in the use of these 
approximations and in the interpretation of results. 

Extensive work has been done using the lowest-order 
model for aircraft flight, vis., the energy model (Refs 1 to 
3 and 18 through 22). In this approximation, the aircraft- 
performance problem is presented in terms of the sum of 
kinetic and potential energies. The control variable is 
either the altitude or airspeed and the state variables are 
the specific energy (kinetic plus potential per unit 
weight), fuel consumed and down- range. A major assumption 
employed v;ith the energy model is that the aerodynamic drag 
can be appro.ximated by its level-flight value. Several 
aircraft trajectory-optimization problems have been handled 
using this model, for example; minimum-time climb, minimum- 
fuel climb, minimum-time turns in a plane, and cruise 
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trajectories have been reported (Refs 15 through 27). The 
relationship between the optimal paths emerging from the 
energy model have been compared with optimal trajectories 
according to more complex models to a certain extent (Refs 
28 through 32). 

The energy model began as an ad hoc approximation (Refs 
1,2,3) and found a theoretical basis in the singular- 
perturbation theory of ordinary differential equations (Refs 
33 through 38). The concept of singular perturbation 
technique in problems of flight mechanics was introduced by 
Kelley (Refs 36, 37 and 38) and applied in different 
formulations by Caiise (Ref 39), Ardema (Ref 32) and others 
(Refs 40 through 45). It has been demonstrated that 
singular-perturbation methods are useful for extending 
energy- state modeling approachs to more general problem 
formulations. These methods constitute a reduced-order- 
analysis approach wherein the system dynamics are separated 
into slow and fast modes. This permits the solution of 
higher-order problems to be approximated in terms of the 
solutions of a series of low-order problems. The singular- 
perturbation procedure has been applied to several problems 
of flight mechanics and has resulted in guidance laws that 
have a nonlinear-feedback form (Refs 46 through 57). These 
schemes may be useful for on-board mechanization. If the 
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present trends continue, this approach appears to hold 
considerable promise . Recently (Ref 58), a new set of state 
variables which offer axitractive replacements for altitude 
and airspeed in singular-perturbation procedures has been 
suggested with a view to enhancing the fidelity of the 
zeroth-order solution as far as possible. Even if the 
results from reduced-order modeling are not directly used 
for on-board implementation, the boundary-layer structure 
and hierarchical ideas from singular perturbations sometimes 
suggest the synthesis of near-optimal guidance schemes (Refs 
59 and 60) . 

Using the assumption that aerodynamic drag can be 
approximated by its level-flight value, a model intermediate 
in complexity between point-mass and energy models can be 
obtained for aircraft in symmetric flight. In this model, 
the flight-path angle is relegated to the status of a 
control variable and the state variables are altitude, 
airspeed, fuel consumed and down-range. An apparent 
advantage of the intermediate vehicle model over the energy 
model is that it can generate possibly realistic path angles 
along optimal trajectories. There is obviously trouble 
ahead with this modeling should the time derivative of path 
angle turn out to be large in optimised maneuvering or, 
worse yet, should the path angle exhibit jump behavior. 
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Although there is some documented research work using this 
model (Refs 61 through 65), it has been less popular than 
the eu'^.'rgy model. Notably, Kelley (Ref 61) used this model 
to analyze a class of time-fuel-range problems. In Ref 62, a 
version of this :..odel has been used to examine fuel-optimal 
paths for transport aircraft. Speyer (Ref 66) investigated 
the solution for aircraft cruise given in Ref 63 using 
higher-order necessary conditions and concluded that, in 
general, steady-state cruise is nonoptimal . He further 
suggested that the optimal-cruise trajectory may be 
oscillatory, and this provided the impetus for research 
efforts reported in Refs 67 through 69. Since oscillatory 
trajectories surfaced in Ref 61 also, in connection with the 

r 

range-optimal climb problem, it was conjectured that these 
might be somehow related. The motivation for the present 
work arose out of these issues. It is of interest to note 
that in Ref 63, a variant of the intermediate vehicle model 
was used to analyze the aircraft-cruise problem. The key 
difference between the model of Refs 61, 65, and Ref 63 is 
that the SinY and CosS' terms arising in the system have been 
replaced by their small-angle approximations. 

With this background, it is the objective of this study 
to investigate optimal time-fuel-range trajectories in 
symmetric flight for aircraft using the intermediate vehicle 
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model. It is stressed at the outset that the investigation 
is oriented in a research spirit of exposing any quirks in 
the model rather than towards any specific application. The 
analysis is based in part upon an exploration of Euler 
solutions for the path-angle-as-control model carried out in 
Ref 61. It happens that under a constant-density-atmosphere 
assumption, the Euler equation for this model can be solved 
in closed form, permitting extensive analytical 
investigation of the solution family. Specifically, optimal- 
range climb and glide problems, and the climb-dash intercept 
problem similar to that in Ref 13 are studied in detail. 
Through second-order necessary conditions it will be shown 
that the minimum-range-to-climb problem (the "steepest 
climb" of Ref 5) has no proper minimum or even a lower 
bound. This problem is found to be ill-posed with the 
vehicle model under consideration, with straight-vertically- 
upward maneuver sequences being furnished by a family of 
paths alternating between upward and downward vertical 
flight and including a limiting chattering member. 
Investigation of the maximum-range-glide problem using 
higher-order necessary conditions is also discussed. For the 
climb-dash inrercept problem, the choice of Lagrange 
multiplier ratios and the selection of an optimal trajectory 
from the Euler-solution family employing the conjugate-point 
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test is given in some detail, An interesting result (Ref 70), 
on the relation between a class of singular optimal-control 
problems and the identically non-regular problems in the 
Calculus of Variations, obtained while analyzing the higher- 
order necessary conditions for minimum time climb and fixed- 
throttle minimum-fuel climb trajectories is presented in an 
appendix. Numerical investigations, including a numerical 
conjugate point test, carried out to verify the conclusions 
from simplified analysis using typical high-performance 
aircraft ( F-4 ) data are also discussed in detail. 

The Euler equation for the time-fuel-range optimal- 
control problem using the intermediate vehicle model appears 
in a nonlinear-feedback form and a cursory analysis reveals 

t 

a scheme for possible on-board implementation; however, 
aerodynamic-drag modeling is questionable. 

1.1 EQUATIONS OF MOTION FOR AIRCRAFT FLIGHT 
The point-mass dynamical model of aircraft flight in three 
dimensions incorporating the assumptions of thrust-along- 
path, zero side-force and flight over flat earth is given by 

V = g[[(T-D)/W| - Siny) (1.1) 

h = V Siny (1.2) 


= (g/v) [ [L cos^/w) - Cosy] 


(1-3) 
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X = g 

L Sin0/(W V Cosy) 

(1.4) 

X = V 

Cosy Cosx 

(1.5) 

y = V 

Cosy Sinx 

(1.6) 

• 

W = Q 


(1.7) 


Additional assumptions embodied in these equations are that 
the throttle is fixed, the variation in weight due to fuel 
expenditure is ignored and winds aloft are zero. The control 
variables in this model are the bank angle 0 and the lift 
coefficient Cj^ or the angle of attack a . The drag is 
modeled as a parabolic function of the lift coefficient; 

""d = ^Do ^ 

The zero-lift-drag coefficient and the induced-drag 

coefficient K are functions of Mach number 

(M) and K = K(M) 

The thrust is a function of Mach number and altitude. 


T = T(M,h) 
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1.2 SYMMETRIC FLIGHT 

In the present work attention will be focussed on aircraft 
maneuvers in the vertical plane. Consequently, the bank 
angle <t> and the flighr path heading angle x can be set to 


zero, resulting in 

V = g[l(T-D)/Wl - Siny] (1.9) 
h = V Siny (1.10) 
y = (g/v) [[L/wj - Cosy] (1.11) 
X = V Cosy (1.12) 
W = Q (1.13) 


Note that the introduction of symmetric-flight assumptions 
has decreased the number of state variables by two. The only 
control variable in this modal is the lift coefficient Cj^ or 
the angle of attack a. 

1.3 INTERMEDIATE VEHICLE MODEL 

In the spirit of model order-reduction, the sweeping 
assumption that aerodynamic drag can be approximated by its 
level-flight value is next invoked. This allows the 
treatment of drag as a function of airspeed and altitude 
only. As a consequence, the flight-path angle y is relegated 
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to the status of a control variable and the equation (1.11) 
can be deleted from the system. The lift coefficient, C^, or 
the angle-of-attack a , previously a control variable, is 
correspondingly assumed to be such as to satisfy equation 
(1.11). An additional order-reduction would lead to the 
energy-state model with altitude or airspeed as the control 
variable. A discussion on the validity of the approximation 
made with regards to the aerodynamic drag is given in Ref 
71. Thus the intermediate vehicle model for aircraft flight 
is given by: 


g[ {(T-D)/W] - Sinr] 

(1.14) 

V Siny 

(1.15) 

V Cosy 

(1.16) 

Q 

(1.17) 


The drag coefficient is now in the form 
Cd = + K 4 W^/(p S)^ 

Note that the Sin!f and Cos2r terms occurring in the system 
equations have been left unmodified although consistency 
demands their replacement by small-angle approximations in 
accordance with the assumption made regarding lift. This is 
done to avoid imposing an artificial bound on path-angle as 
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will be required otherwise to obtain sensible results along 
certain optimized paths IRef 17). It should be noted that 
the approximation carried out here cannot properly be 
classed as a singular perturbations procedure as in Ref 38. 
It must, therefore, be accepted at the outset that those 
solutions which violate the assumption that i is small along 
substantial portions of the optimized trajectory should be 
interpreted with caution. In the model given above, the 
throttle has been assumed fixed. It will be seen in 
subsequent portions of this work that the throttle control 
can be reinstated without excessive difficulty. However, 
this tends to complicate the analysis somewhat and is 
treated in the present work only for completeness. 

1.4 j OPTIMAL -CONTROL PROBLEM 

To state the optimal-control problem concisely, it is 

desired to find the control history ^(t) such that the 

sysrem described by (1.14) through (1.17) be transferred 

from initial state (V ,h^,j: ,W ,t ) to the final state 

o o o o o 

( h,-, x^,W^, t^) , not all of whose components are specified, 
such that a functional P(V_,h,r/X^, W,, t,.) takes on a minimum 
value . 

From a practical viewpoint, the time-range optimal- 
control problems are of main interest since minimum-fuel 



13 


problems with fixed throttle rarely occur in applications. 
Therefore, major attention is focussed on the analysis of 
full-throttle, minimum-range climb (Steepest climb of 
Ref, 5), power-off maximum- range glide and the climb-dash 
intercept problems. Out of these three, the climb-dash 
intercept problem requires some explanation. This trajectory 
occurs in aircraft pursuit-evasion problems (Ref 13, 38) as 
a transient leading to the dash point or the high speed 
point on bhe level flight envelope. In general, the time 
spent during the climb is much smaller than the ti;ae spent 
in the dash. This maneuver normally ends with a terminal 
transient as mentioned in Ref. 13. The treatment in this 
work, however, will not include any discussion of this 
terminal transient. 



Chapter II 

OPTIMAL CONTROL WITH INTERMEDIATE VEHICLE MODEL 

Introduction of the assumption that the aerodynamic drag can 
be approximated by its level-flight value in the equations 
of motion for symmetric flight leads to a model intermediate 
in complexity between point-mass and energy models. The 
chief advantage gained in employing this approximation is 
the reduction in model order offering analytical and 
numerical simplicity. The control variable in this model is 
the flight path angle Jf and the state variables are the 
airspeed, altitude, down-range and fuel consumed. 
Presumably, the solutions emerging from this modeling may 
have validity at least along central portions of the 
optimized path. It should be noted, however, that the 
intermediate vehicle model permits jumps in the path-angle 
and, hence, the assumption of small Z is violated along the 
optimized trajectory. Consequently, care should be 
exercised in the interpretation of the results. 

The optimal-control problem under consideration is the 
minimization of a function of the state-variables and time 
final values, subject to the differential constraints given 
by the equations (1,14) through (1.17), with satisfaction of 
any boundary conditions that may be imposed. To this end, 
the variational Hamiltonian function is formed as 


14 
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H g [{(T-D)/Wj - Siny] + V Siny 

+ Xj^ V Cosy + x^ Q (2.1) 

The Euler-Lagrange equations are 

V " "^V 

- x^ Cosy - X^ 32. <2.2) 

Xh = -X ^i-<T-D) - x~ ^ (2.3) 

X^ = 0 (2.4) 

X^ = 0 (2.5) 

The optimality condition 3H = 0, is given by 

ay 

-X^ g Cosy + X^ V Cosy - X^ V Siny =0 (2.6) 

In these equations, Siny an"’ Cosy terms have been retained 
unmodified, although it is evident that consistency demands 
their replacement by small~angle approximations in 

accordance with the assumption made concerning the 

aerodynamic drag. This has important bearing on the question 
of imposed boundary conditions, as will be seen 

subsequently. The analysis presented in this section closely 
follows that of Ref 61. 
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Equations (1.14) through (1.17) and the equations (2.2) 
through (2.6) together with specified boundary conditions 
and appropriate transversality conditions form the two- 
point-boundary-value problem, which needs to be solved to 
obtain the optimal control t{t) . If desired, a numerical 
solution can be obtained for this system. Be that as it 
may, the approach adopted here is in a different direction. 
It can be noted from the Euler-Lagrange equations (2.2) 
through (2.6) that the only costates that are variable in 
the Hamilronian are those associated with the airspeed and 
altitude, vit., Xy and X^^ . If these can be eliminated 
using algebraic manipulations, then a first integral of the 
motion, i.e.,the variational Hamiltonian, can be obtained in 
a closed form. 

In the following, the time-derivatives of the 

optimality condition (equation 2.6) will be used to 

eliminate the costates X,, and X. in favor of the control 

V n 

variable y and its derivatives. Note that this procedure is 
somewhat formal since the derivatives of the path angle may 
not exist at certain points on the optimized trajectory. 
Using equation (2.6) 


(2.7) 



17 


ORiGf|\’Al r/ 

OF POOR QUALnV 


Differentiating the above expression once with respect to 
time, 

Kj = V (X, - X Tanir) + V (X, - X Z Sec^Z) (2.8) 

V— 'hx — 'hx ' 

<3 g 

Now, substituting for X^ in equation (2.3) from equation 
(2.7) 

K = - 1 (K - Tanar) MT-D) - X- ag (2.9) 

^ w ^ ^ ^ aK 

Similarly, the equation (2.2) becomes, 


K; - - 1 (K - Tany) i_(T-D) - X, Siny 
^ w ^ ^ 3 V ^ 

- X Cosy - x^ ^ (2.10) 

^ av 

Substituting next for X^ from equation (2.8) in the equation 
(2.10), one obtains a second expression for Xj^ as 


= - ^h ^ 

) + 1_(T-D)1 

W L V 

av J 

- K,\ 9 ‘ y 

Sec^y - g Tany 

^ 1 V Cosy 

W 


i_(T-D) ll- X- g a£_ 
av V av 


( 2 . 11 ) 


t,quations (2.9) and (2 11) can now be used to eliminate X^^, 
resulting in the following expression relating Xj^, X^, 


and 
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X., f8_[V(T-D)) - ^ i_lV(T-D)l I 
W [ah V 3V J 

- X fTany f ^[V(T-D) ] - £ ^[V(T-D) ] 
w^'L \ ah V av 

+ g - y sec^y] + = 0 

V cosy J Idh V avj 


( 2 . 12 ) 


The expressions (2.7) and (2.12) for X^ and Xj^ is next 
substituted in the variational Hamiltonian (2.1) with the 
following expression resulting 


Cosy H 
- Cosy x^ 



(/a - g L_^ [V(T-D) ]) 

V av/ I 

Q* i - 5. IV(T-D)/Q] 
) V av/ 


- 2 . 

V 


av/ 


[T-D) 


(T-D)V y 

Cosy 


0 


(2.13) 


It may be noted that 


/a - g 

i-) 

[ 

) = L. [ ] 


lah V 

av/ 


ah 

E = Constant 


where E = h + , the specific energy. 

2g 

The first integral of motion (2.13) is now independent 
of the multipliers X^ and If the required data is 


available, this expression together with the differential 
equations (1.14) through (1.17) can be numerically 
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integrated to obtain optimal trajectories} provided the 
initial value of path-angle is known. Such solutions are 
given elsewhere in this work. Rather than attempting to 
generate numerical solutions immediately, it is desirable at 
this point to examine specific cases to gain better insight 
into the nature of trajectories emerging from the expression 
(2.13) . 

For completeness, the inclusion of throttle control in 
the above formulation is next considered. One of the ways to 
do this is to introduce a multiplicative parameter, n , on 
thrust, T, and fuel-consumption rate, Q, in the stace-Euler 
system, with 


0 < T1 < 1 

The thrust and fuel-consumption rates appearing in these 
equations are then their maximum values at each altitude and 
airspeed. In such a situation, the Minimum Principle (Ref 
24) requires that the throttle setting be chosen at each 
instant such that. 

n = 0, if 9 .T+X~qJ >0 

n = 1, if ^X^ 

The multiplier X^ required in the switching function can be 
determined from equations (2.7) and (2.12). Note that a 
singular arc arises whenever 
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-2 ^ 2 - ° 

w 

over a nonzero interval of time. In this case, the control, 
T), is determined by equating successive time derivatives of 
the switching function to zero. It is also known that the 
Generalized Legendre-Clebsch necessary condition, sometimes 
known as the Kelley-Contensou test (Ref 72), should hold 
along the minimizing singular arcs. This will not be 
pursued any further in the present work and, in the 
treatment that follows, the throttle will be assumed fixed 
at its maximum value. 

2.1 PROPERTIES ^ lliE SOLUTION 

In order to investigate the implications of expression 
(2.13), three distinct cases will be considered in the 
following. These are minimum-time-to-climb, minimum-fuel-to- 
climb and range-optimal climb and glide trajectories. To 
illustrate the approach, consider, for example, that the 
final values of range and time are open. Then the 
transversality conditions X = 0, H = 0 apply and the 
optimization problem is a trade-off between the final values 
of fuel-consumed, altitude and airspeed, the maximum or 
minimum value of one of these variables or some function of 
these variables being sought without regard to the final 
values of range and time. 
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2.1.1 Minimum-Time-To-Climb Trajectory 

If time to climb is to be minimized, then the Hamiltpnian H 
has a value of -1 ( Ref 23,61 ), = 0, = 0 and the 

expression (2.13) reduces to 

Cosy ( /a -ga \ [V(t-d)1 I = o (2.14) 

i 1 3h V avy ) 

which may also be written as 


Cosy ^ [V(T-D) 
ah 


= 0 


(2.15) 


E = Constant 

The equation (2.14) can be satisfied either by Cosy = 0, 
vertical flight, or by the vanishing of the expression in 
brackets. This expression is the derivative of the so-called 
"excess power" in the flight performance literature with 
respect to altitude with the specific-energy held constant. 
The minimum-time-to-climb trajectory may be conceived as a 
path in the altitude-airspeed chart which passes over each 
specific-energy curve at a local maximum of specific excess 
power along that curve. This trajectory carries the name 
"energy climb" (Ref 1) due to its close association with the 
specific-energy concept. It is appropriate to point out that 
the vertical-flight result would not have been obtained if 
the Siny and Cosy terms in the system had been replaced by 
small-angle approximations. 
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Thus the solution to this, or any, altitude-airspeed- 
time optimal problem is made up of vertical climbs, vertica] 
dives and energy climbs pieced together in proper sequence; 
a result noted by several researchers (Refs 1 through 6 and 
15, IS, 19 and 61) . 

2.1.2 Minimum- Fuel-To-Climb Trajectory 

To generate a minimum-fuel climb path, with final values of 
time and range open, the natural boundary condition X = 0 
and H = 0 can be used in equation (2.13), with the following 
result 

Cosyj [V(T-D)/Q] 5 = 0 (2.16) 
X I 3h V av / j 

Just as in the energy-climb case, the expression (2,16) can 
be satisfied either by Cosy = 0, vertical flight, or by 
vanishing of the bracketed expression. This quantity is the 
derivative of the excess power per unit fuel-flow with 
respect to altitude, with specific energy held constant. 
Hence, the solution of minimum-fuel-to-climb is made up of 
vertical climbs, vertical dives and a path in the altitude- 
airspeed chart which passes over each specific energy curve 
at a local ma.ximum of the excess power per unir fuel-flow 
along that curve. This trajectory is sometimes referred to 
as an "economy climb". 
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2.1.3 Range-Optimal Climb and Glide Paths 
If the range is to be minimized or maximized with final time 
and fuel unspecified, then = ± 1 and H = X^ = 0, and a 
first-order differential equation emerges from expression 
(2.13) for path inclination as follows 

/ 3 - g 3 Vt-D) - (T~D)V Y = 0 

Vsh V av/ Cost 


(2.17) 

Note that the Z term would not have appeared in the above if 
small-angle approximation had been used in the equations of 
motion. This has a significant impact on the nature of 
optimal trajectories for this problem, as will be seen later 
in this chapter. The system consisting of equations (1.14) 
through (1.17) and (2.17) generates a trajectory family for 
the range-optimal problem. The possibility of obtaining 
analytical solution to this system for the case of thrust 
and drag as arbitrary functions of altitude and airspeed is 
remote. As a result, the interpretation of expression (2.17) 
is not as straightforward as it was in the minimum-time-to- 
ciimb and minimum-fuel-to-climb cases. 

However, introducing the assumption of constant-density 
atmosphere with thrust and drag dependent on airspeed only, 
it is feasible to obtain an analytical solution to this 
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system ( Ref 61 ) . This assumption has been employed in the 
past for the analysis of some optimal aircraft trajectories, 
see Refs 73 and 74, for example. Results obtained from this 
analysis may perhaps be characteristic of more general 
solutions as well. 

The equation (2.17) can be rewritten as 


Since time is not of particular interest in the problem, it 
may be eliminated in favor of airspeed as an independent 
variable. In the following/ several transformations of the 
independent variable are carried out without attention to 
monotonicity requirements, the thought being to fit solution 
segments obtained into families in due course. The 
temptation of range as independent variable will be avoided, 
however, in anticipation of purely-vertical-motion segments. 
In the interest of brevity, designate the acceleration 
variable \i = {T-D)/W, 

® ^ ^ 

^ V 3Vi 

( 2 . 19 ) 


(]i - Siny) = _V 

gy 


= V ) ( L. ~ ^ lA 

Cosy (T-D) ly V 
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With altitude dependence of thrust and drag supressed, the 
path angle Y is determined as the solution of the first- 
order differential equation 

1 dy (SinJf - ij) “ A. ^ (2.20) 

Cosy dV y dV 

Further simplification is obtained by another change of 
independent variable, this time from airspeed V to y 

1 ^ (Siny - y) = _1, (2.21) 

Cosl? ay y 

If the r61es of independent and dependent variables are now 
regarded as reversed, this equation takes the form 

dy + 1 y^ - y Siny = 0 (2.22) 

dy Cosy Cosy 

which is a form of the Bernoulli differential equation 

y* + f,(y) y^ = 0 (2.23) 

dy ^ 

with 3 = 1. According to Kamke, ( Ref 75 ), this equation 
has the solution 

_1_ = E(y) dy (2.24) 

y J E{H) 

where 

/fo(y) dy 

E(y) = e (2.25) 

with the identification of f^ and £2 as 
= 


_1 

Cosy 


(2.26) 
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f2(ir) = 


/ 


SinJf dy 


Cosy 


The solution ( 2 . 24 ) becomes as follows 


In Cosy 

E{y) = e = Cosy 



1 

V = siny + c Cosy 


(2.27) 


(2.28) 


(2.29) 


Before expressing this relationship in the form y = 

r 

y(v), one may relate the integration constant C to 

equilibrium values of p and y corresponding to unaccelerated 
flight. Such values may be designated with a superscribed 
bar : 

il = Siny (2.30) 

C = Coty (2.31) 


The solution can then be expressed as : 


Sin7 Siny + Cosy Cosy - — 

P 

or as 



(2.32) 


(2 .33) 
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Here ]! is the value of v in unaccelerated flight and 

? = Sin'^it (2.34) 

In Fig. 1, the solution (2.33) for the range-optimal 
problem is illustrated for various values of u. Note that 
the range of the flight-path angle Y has been restricted to 
il80° in this plot, 

The solution family in Fig. 1 is not related 

sufficiently closely to the state variables in the system to 
be very illuminating. This may be remedied in part by 
reinstating airspeed as the independent variable. To this 
end, consider first the case of unpowered flight with a 
quadratic distribution of power variable vs airspeed as 
shown in the upper part of Fig. 2. The solution family 

generated by substituting the assumed quadratic in equation 
(2.33) for various values of parameter y is sketched in the 
lower portion of the figure. The point identified by a 
circle corresponds to steady flight at minimum-drag airspeed 
or glide. The arrows on the solution family indicate 

the direction of increasing time, determined from the 
differential equation for airspeed, viz.. 


V 


g( y - SinY ] 


(2.35) 
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The family of transients leading to and from the L/D 

max 

point can also be seen in this figure. Other paths to the 
right and left of this plot appear to be connected with 
transient dive and zoom solutions which do not involve 
prolonged operation at or near the minimum-drag airspeed. 
Thus, according to the present modeling, the solution to 
maximum-range glide is comprised of the following. A 
transient path leading from the initial condition to the 
L/D flight condition, flight at the minimum-drag speed, 

InaX 

and another transient terminating at the final boundary 
conditions; all this provided that the terminal conditions 
are attainable in gliding fxight. If the initial condition 
does not lie on a transient path leading to the L/D 

max 

point, a jump in the path angle JT is indicated, which places 
the vehicle on the appropriate trajectory. Similarly, if the 
terminal boundary condition is off a transient path leading 
away from the minimum-drag airspeed, another jump in path 
angle is required. Therefore, with intermediate vehicle 
modeling, the maximum-range glide path, in general, will 
contain points at which the path-angle jumps. Consequently, 
the assumptions made with regards to the magnitude of the 
derivative of the path angle will be violated along the 
optimized path unless the specified boundary conditions 
happen to be on a transient path and sufficiently close to 
the equilibrium flight condition. 
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When a positive margin of thrust over drag exists, a 

family of oscillatory solutions corresponding to minimum- 

range climb can be generated, again by assuming a quadratic 

relation between the acceleration variable and airspeed, as 

depicted in Fig. 3. The innermost circled point in this 

chart corresponds to flight along (T-D) . A family of 

max 

oscillatory solutions about this point appears with an 
outermost limiting member along which the path angle 
switches between ± 90°. Interpretation of these 
trajectories, in the light of results given in Ref 5 for uhe 
"steepest climb" problem, is not straightforward. However, 
an examination of second-order necessary conditions to be 
presented in the next chapter will reveal that the "steepest 
climb" problem does not possess a minimum or even a lower 
bound, when a positive margin of thrust over drag exists. 

With the experience gained from the analysis of various 
special cases, the next obvious step is to tackle the more 
general time-fuel-range problem. If one chooses = -1, X^ 
= 0, and a fixed value of H (to be determined), the 
expression (2.13) is the Euler equation for the problem of 
maximising range with a fixed final time. Similarly, if X = 
-1, H = 0 and X^ is fixed, the Euler equation for maximum- 
range trajectory with fixed final value of fuel is obtained. 
It may be noted that the range-maximization problem without 
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time or fuel • -istraints does not have a proper maximum or 
an upper bound. It is clear that fixed- throttle range-fuel 
optimum problems are of little interest in practical 
situations. Hence, in the next section, attention will be 
focussed on the problem of maximizing the range with a 
.specified final time. 

2.2 T IME-RAMGE OPTIMAL PATHS 

With = -1, X^ = 0, and a fixed value of H, expression 
(.2.13) is the Euler equation for the time-range problem, In 
its original form, this equation is not amenable to detailed 
analysis due to its nonlinear nature and the appearance of 
arbitrary thrust and drag functions. However, it happens 

f 

that, if one employs the constant-density-atmosphere 
assumption, analytical solution to the general expression 
(2.13) can- be obtained using the variatioh-of-parameters 
technique (Refs 76 and 77) on the solution for the range 
problem. Required transversality conditions can then be 
imposed to study speci'fic cases. 

Under the constant-density-atmosphere assumption, 
equation (2.13) can be written as 

- g( cosy ^ - cos^'y (h/x ) (i/v^v) ±_{Vv) 

\ u dV dV 

+ Cos^y (X^/Xj,) (QVV^li) d_[Vw/Q]l 


(2.36) 
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The independent variable is nov; changed from time to 
airspeed, resulting in 


dy (SinJT - y) = CosYdu 
dV M dV 

- Cos^y (l/v"y) d_(Vy) 

^ dV 

+ Cos"y (QVV*y) d_[Vy/Q] 

” ^ dV 

(2.37) 


Rearranging, 


dy (Siny - y) ■' Cosy dy = - Cos*y T (H/X ) (1/V“y) d (Vy) 
dV y dV L ^ dV 

“ (QVV^y) d [Vy/Q]“| 

dV J 


(2.38) 

The analytical solution for the first-order differential 
equation (2.38) with the right members set to zero is the 
expression (2.29). The equation (2.29) may now be 
differentiated with respect to airspeed, 

- i_ ^ = (Cosy - c Siny) ^ ^ Cosy 

dV dV dV 


(2.39) 



32 


ORIGINAL PAGE 
OF POOR QUALITY 


Note that C is no longer a constant, but a function of the 
independent variable V. Substituting for y in equation 

(2.39) from (2.29), 

.1 ^ ■ f(cosy - c siny) ^ ^ Cosy 1 

]i dV L dV dV J 

(Siny + C Cosy) 


(2.40) 


Employing the expression (2.40) in (2.^38) and carrying out 
algebraic manipulations, one obtains a first order 
differential equation for C as 

Cos^y ^ = - Cos^y (H/X ) (1/V*y) d_(Vy) 

Siny + C Cosf dV ^ ' dV 

+ Cos^y (X~/X ) (QVv*y) ±_m/Q] 

fV JX7 


(2.41) 


Simplifying (2.41) using equation (2.29) 

-(H/X ) I (1/V*y) + (1/Vi/) 

dV ^ 1 dV I 

+ ! (Q/V“Vi) + (Q/Vy^) ^ - (1/Vy) ^1 

^ 1 dV dV / 


(2.42) 


The Quantities within the brackets can be identified as 


- ^[1/Vy] = (1/V*y) + (1/VyM ^ 
dV dV 


(2.43) 
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and 


- ^IQ/Vv] = > ( 1 /Vv) ^ + (Q/VvM ^ + (Q/V*y) 
dV Hv dV 


(2.44) 


from which 


^ = (H/X^) ^[1/Vy] - (Xjv/X ) ±_[Q/V]i] 
dV ^ dV ” ^ dV 

Equation (2.45) is readily integrated to yield 
C = H. ^ + C. 


X Vv 

X 


-W - 

x” vy 


(2.45) 


(2.46) 


Where is an arbitrary constant. 

Hence, the solution fox. time-range- fuel optimal control 
problem with altitude dependence of thrust and drag 
suppressed is 


1 = Siny 

y 


+ r H 1 

L X, Vy 


X ■* 


Cosy 


(2.47) 


To express the above result in the form y = y(y), the 
integration constant: needs to be related to equilibrium 

values of y and H corresponding to unaccelerated flight. 
Unlike the siruation in the range optimal problem, the 
interpretation of the solution (2.47) is not 
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straightforward, partly due to the involvement of time-range 
and fuel-range multiplier ratios in this expression. 

As discussed earlier, fixed-throttle minimum-fuel 
trajectories are of limited interest in practical 
situations. Hence, in the following, primary attention will 
be focussed on the time-range optimal-control problem. This 
permits the deletion of a term in (2.47) resulting in 


1 = Siny 

y 


+ H 1 + 

T vT 1 

m ^ • 


Cosy 


(2.48) 


To investigate the range of permissible H/X values, the 
equilibrium flight conditions obtained from equations (2.13) 


and (2.35) are next examined. Setting V and y equal to zero 
in these expressions results in 


Siny = u (2.49) 

Cosy ^(Vu)l (2.50) 

dV J 

Combining equations (2.49) and (2.50), an equation for the 
multiplier ratio H/X corresponding to equilibrium flight 

conditions, in terms of airspeed and the acceleration 

variable emerges as 

H. = V* ^ yi-y^" ^(Vy) 

X dV L dV 

The equation (2.51) may be evaluated in the range of 




(2.51) 


H = dw / 
T dV 


airspeeds of interest to obtain permissible values of the 
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multiplier ratio H/X„. In Fig. 4 a typical parabolic 
distribution of the acceleration variable y vs airspeed is 
illustrated. Employing this distribution in equation 
(2. Si), one can identify the permissible values of H/X . A 

A 

plot resulting from this calculation is shown in Fig. 5. In 

this figure, three distinct regimes may be identified. H/X 

values to the left of (T-D) point (marked A in Fig. 5) 

max 

are positive, while those between the (T-D) point and the 

max 

V(T-D) _ point (marked B) have a negative sign. All H/X 
IHcLX X 

values to the right of the airspeed for V(T-D) are 

positive. Any of these values may be used to evaluate the 
arbitrary constant as follows. As in (2.30), 

y = Sin? (2.52) 


V = V 

Equilibrium value of E 


Coty = H 1 + C 

Or 


~ 

Coty 

- H 1 



X Vil 

X 

Using (2. 

55 ) in 

(2.48) , 

1 = Siny 

.Th 

j 1 - 1 1 + Coty 1 

y 

X 

L X 

/ Vy W t J 


Cosy 


(2.53) 


(2.54) 


(2.55) 


(2.56) 


Now, put 
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i = H U - 1 + Cot? 

1 Vv vp I 

and using a well known trignometric identity, 

y = Tan'^ Fi-V Cos’^r ^ 1 (2.57) 

L A J L v'a^ + 1 J 

Expression (2.57) is the solution of the Euler equation for 
the time-range optimal control problem with altitude 
dependence of v suppressed. In Fig. 6, \;he solution (2.57) 
evaluated using a representative H/X value from the first 
equilibrium regime is shown. From this figure it may be 
observed that solution family is oscillatory about the 
equilibrium point. A similar family of solutions is 
obtained if the multiplier ratio H/X from the second 
equilibrium regime is used, as illustrated in Fig. 7, Using 
H/X values from the third equilibrium regime, i.e.,, values 

X 

to the right of V(T-D) in Fig. 5, results in the solution 
family given in Fig. 8. In this figure a paix" of 
transients leading to and away from the equilibrium point 
may be observed. Other transient paths appear to be 
connected with those flight conditions which do not involve 
prolonged operation at or near the equilibrium point. It is 
of interest at this juncture to point out the similarity of 
Figs. 6 and 7 to Fig. 3. The solutions in Fig. 8 bear some 
resemblance to the plot of Fig. 2, 'It will be seen 
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subsequently that the oscillatory paths of Fig. 6 and 7 
generally fail the conjugate-point test for durations 
exceeding one period and furnish a minimum only for short- 
duration problems. According to Fig. 8, then, "long" 
optimal time-range trajectories are made up of a transient 
path leading to the equilibrium conditions, flight at the 
equilibrium condition and a transient path leading to the 
terminal boundary conditions. If the boundary conditions are 
off the transient paths, jumps in path-angle take place. 
Hence, it is clear that the assumption made with regard to 
the magnitude of f will be violated on segments of the 
optimal time-range trajectory unless the boundary conditions 
happens to be on a transient path and sufficiently close to 
the equilibrium flight condition. 

2.3 CONCLUSIONS 

In this chapter a study of the time-range-fuel optimal 
control problem for aircraft using the intermediate vehicle 
model was presented. Employing the constant-density- 
atmosphere assumption, the range optimum climb and glide, 
and the time-range optimum problems were studied. One notes 
that the range problem has oscillatory solutions when a 
positive margin of thrust over drag exists. With zero 
thrust, the solution obtained is the flattest glide with a 



38 


family of transients leading to and from the L/D point. 
For the time-range problem, values of the multiplier ratio 
H/X„ to the left ( low-speed side ) of the V(T-D) point 
produce oscillatory solutions. With H/X values chosen to 

X 

the right of airspeed corresponding to a solution 
family consisting of a set of transients leading to and from 
the equilibrium point, defined by the particular choice of 
the the multiplier ratio, is obtained. Optimality of these 
Euler solutions are investigated in the next chapter using 
second-order nece.ssary conditions. To verify the conclusions 
arrived at in this chapter, numerical computations are 
carried out using the data for the F-4 aircraft in Chapter 
IV. 



Chapter III 

SECOND - ORDER NECESSARY CONDITIONS 

It is known that the sufficient conditions for a functional 

J(x) to have a "weak" minimum for x = x are that the first 

* 

variation 6J(x) vanishes for x = x and that the second 

, * 
variation 6*J(x) be strongly positive for x = x (Ref 78). 

Investigation of the conditions for vanishing of the first 

variation leads to the Euler-Lagrange equations and 

transversality conditions. The positivity requirement of 

second variation of the functional J(x) leads to the 

Legendre-Clebsch and Jacobi necessary conditions. If both 

the strengthened Legendre-Clebsch condition and strengthened 

Jacobi test are satisfied by an extremal, it can be shown 

that the second variation of the performance index J(x) is 

positive for sufficiently small variations (Refs 78 through 

86). Thus, the functional J(x) has a weak relative minimum 

if the extremal satisfies 1) the transversality conditions 

2 ) the strengthened Legendre-Clebsch necessary condition 

and, 3) the strengthened Jacobi condition. In addition, if 

the e.xtremal satisfies the Weierstrass condition, it affords 

a strong relative minimum for J(x). In optimal-control 

problems, when the control vector is not subject to 

inequality constraints, it can be shown that the Pontryagin 
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maximum principle is equivalent to the Weierstrass condition 
(Refs 24,79,86 and 87). Hence, in these problems if the 
extremal satisfies the transversality conditions and the 
strengthened Legendre-Clebsch necessary condition together 
with the strengthened Jacobi condition, it provides a strong 
relative minimum for the performance index J{x). For a 
relative maximum the preceding sufficiency test applies if 
the inequality associated with the strengthened Legendre- 
Clebsch test is reversed. 

In the following, families of Euler solutions obtained 
in Chapter II for range-optimal climb and glide, and time- 
range optimal problems are examined using the Legendre- 
Clebsch and Jacobi necessary conditions. It will be seen 
that the constant-density-atmosphere assumption enables the 
analytical treatment of the Jacobi test for these problems. 

3.1 LEGENDRE-CLEBSCH NECESSARY CONDITIONS 

From Euler-Lagrange equation (2.6), with the fuel multiplier 
set to zero, 

31^ = - g Cosy + V Cosy “ V Siny (3.1> 

3y 

and 


= (X^ g - X^ V) Siny - ^ Cosy 


ay^ 


(3.2) 
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Setting the Left member of equation (3.1) to zero as 
required for a stationary minimum of the Hamiltonian leads 
to 


Tany = (Xj^ V - Xy g)/X^ V 


Or 


Tany = (X^ g - V)/(- X^ V) 

From (3.3), then 

• Siny = (X^ V - Xy g) 0 

\/(X, V - X,, g)^ + X* ^ 
' h V ' X 


(3.3) 


(3.4) 


and 


Cosy = X^ V 0 (3.5) 

^ 

v/(\ V - x^ g)* + X*^ V* 

f 

where o = ±1. Employing (3.4) and (3.5) in (3.2), the 

sign of o can be determined. 

Next, the transversality conditions for the range 
problem, viz., X = 1 for range minimization and X = -1 for 

X X 

range maximization is used in (3.2), with the following 
result . 

X^ = 1, liH > 0 if y lies in the 

^ ay^ 


second or third quadrant 


(3.6) 
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X 


X 


- 1 , 


3*H < 0 if y lies in the 

TF’ 

first or fourth quadrant 


(3.7) 


From (3.6) it is clear that, with no restrictions on 
path angle Jf, the minimum-range-climb trajectory is that 
which maximizes the range in the negative direction, a 
result which ’i perhaps obvious. The implication is that 
with no constraint on the final value of time or fuel, the 
•'steepest climb" problem does not possess a minimum or even 
a lower bound. 

Attention is drawn to the solution to this problem 

given by Miele (Ref 5) using the Green's theorem device. 

According to Ref 5, the optimal trajectory for the "steepest 

climb" problem consists of a central path flown along the 

(T-D)^^j^ locus in the airspeed-altitude chart with vertical 

climb/dive transitions at the ends to meet the boundary 

conditions if they are off the (T-D) path. There is an 

max 

important difference in vehicle modeling from that of the 
present work which should bs noted as a key to resolving 
disparities between the character of optimal paths emerging: 
The analysis of Ref 5 in essence replaces Cosy in equation 
(1.16) with unity so that the problem solved is maximum 
altitude in a given distance (arc length) rather than in a 


given range. 
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Consider, next, the imposition of limits on path-angle 
y, say -90° Jf < 90°. In this case, one observes that by 
alternating between vertical-climb and vertical-dive paths, 
the range- to- climb can be made identically zero. This is a 
consequence of the intermediate vehicle modeling in which 
there is no limit to the steepness of climb. Examination of 
energy rate 

E = V(T-D)/W 

shows that it is positive as long as there exists a positive 
margin of thrust over drag. Since energy can be traded-in 
for altitude or airspeed, any combination of these can be 
reached in zero range using a composite vertical climb/dive 

•t 

trajectory. Additionally, there is no preferred location 
within the flight envelope for switching from vertical climb 
to vertical dive and vice versa. It is concievable, 
therefore, that in the limiting case, the minimum-range- 
climb path will tend to a "chattering arc" with infinite 
number of switchings between the climb and dive flight 
segments . 

Chattering arcs can also arise in the minimum-time 
climb and minimum-fuel climb problems. Consider, for 
example, a parabolic distribution of specific excess power 
vs airspeed as shown in Fig. 9. In Chapter II it was shown. 
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for the minimum-time climb problem, that the Euler equation 

(2.14) can be satisfied either by vertical flight, Cosy = 0, 

or by flight at the airspeed corresponding to the local 

maximum of specific excess power (V in Fig. 9) at each 

★ * * 

energy level. Choosing a pair of airspeeds V and V about 
V, a trajectory consisting only of vertical climb/dive 
flight segments switching between V and V can be 

constructed. The choice of these airspeeds is arbitrary, as 
long as the specific excess power V(T-D)/W > 0, since this 

•k kk 

will ensure a net gain in energy. If V and V are chosen 

just to the left and just to the right of V, minimum-time 

performance is approached. In the limit as they approach V, 
a chattering solution is defined along with a lower bound on 
the time index. For general end conditions, straight-up 
and/or straight-down transients precede and follow the 
chattering subarc in a composite solution. 

On the other hand, the Legendre-Clebsch necessary 
condition is met in the strengthened form along the maximum- 
range-glide and time-range optimal trajectories for the 
values of path angle y in the first and fourth quadrants. 

It may be noted at this point that in the cases of time 

or fuel-minimization problems with range open, the Legendre- 

Clebsch necessary condition is met only in weak form along 
central arcs and hence, these trajectories fall in the class 
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of singular extremals. In general, demonstrating sufficient 
conditions for singular extremals to be minimizing is a 
nontrivial task, though extensive results exist (Refs 72,88, 
^ and 89). However, a theorem originally due to Mancill (Ref 
90), and widely applied (Refs 4 through 7 and 91) in 
problems of flight mechanics, enables one to establish the 
sufficient condition for a strong relative minimum, along 
central arcs, for minimum-time-climb and minimum- fuel -climb 
problems as 

(i) Minimum- time-to-climb 

^[V(T-D)]| =0 (3.8) 

3h |E = Constant 


ii [V(T-D)1 I < 0 (3.9) 

3h* |e = Constant 

(ii) Minimum-fuel-climb 


3_[V(T-D)/QJ = 0 (3.10) 

3h E = Constant 


[V(T-D)/Q]| < 0 (3.11) 
3h^ |E = Constant 

In Ref 72, the inequality (3.9) was obtained using the 
Generalized Legendre-Clebsch necessary condition in 
conjunction with energy modeling. A question that arises 
naturally here is the relation between Mancill 's work and 
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the Generalized Legendre-Clebsch necessary condition. In an 
appendix, this is discussed in detail. Suffice it to state 
here that, for problems containing two non-ignorable state 
variables and one control, under a smoothness hypothesis, 
this theorem yields a sufficient condition for a strong 
relative minimum. 

3.2 JACOBI ' S NECESSARY CONDITION 

In the earlier section, it was shown that the Legendre- 
Clebsch necessary condition is met with a margin for 
maximum-range glide and the time-range optimal problems. 
Consequently, the Euler solutions for these problems are 
optimal for initial and terminal boundary conditions 
sufficiently close together. For extremals of finite 
leng+'h, however, the task of ensuring that the second 
variation is non-negative for admissible neighboring paths 
leads to the accessory-minimum problem in the Calculus of 
Variations. This in essence boils down to a search for a 
system of admissible variations, not identically zero, which 
offer the most severe competition in the sense of minimizing 
the second variation. If a system of nonzero variations can 
be found which makes the second variation zero, then it is 
clear that a neighboring path is competitive and that the 
test extremal furnishes at best an improper miniminn and at 



worst a merely stationary value (Ref 92). The first value of 
the independent variable x = x* ^ ^o which such a non 
trivial system can be found defines a conjugate point. One 
approach to the accessory-minimum problem in optimal-control 
theory consists of expanding the augmented' criterion ( the 
variational Hamiltonian) to second order and all constraints 
to first order to obtain a matrix-Riccati equation (Refs 23 
and 87). Existence of a bounded symmetric -matrix solution to 
this Riccati equation, then, determines -the existence of a 
conjugate point. An alternate procedure for the conjugate- 
point test was suggested in Ref 92, involving the 
examination of the rank of a matrix of variations of the 
states with respect to the initial values of the costates. 
This procedure will be adopted in the present work. It needs 
to be stressed that both these approachs lead to the same 
conclusions and choice between them is entirely a matter of 
convenience. 

According to the analysis of Ref 92, for the Mayer 
problem, the rank of the matrix of variations of states and 
the multiplier corresponding to the state being minimized 
with respect to the initial values of the costates, 
evaluated along the test extremal, viz. 
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the rank of 
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provides the criterion for the existence of a conjugate 
point. If the rank of the matrix (3.12) drops at any point 
along the test extremal, it is indicative of the existence 
of a conjugate point. 

It is evident that numerical procedures may be set up 
to evaluate the elements of the Jacobi -test matrix, see Ref 
93 for example. Details of one such scheme will be 
discussed in the next chapter. The intention is to test the 
Euler solutions obtained for maximum- range glide and time- 
range optimal problems with altitude dependence of thrust 
and drag supressed, for conjugate points. In view of the 
particularly simple form of the conjugate-point test for 
these problems, it seems reasonable to attempt to obtain 
analytical approximations for the partial derivatives 
required in the Jacobi test matrix. 
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OF POOR QUALITY 


3.2.1 Conjugate-Point Test for MaKimum-Range Glide 
The maximum- range glide problem may be thought of as the 
maximization of the final altitude with zero thrust, for a 
specified range with time and fuel open. Since it is 
convenient to work with range as the independent variable, 
the equations of motion (1.14) through (1.16) are rewritten 
as 


V’ = g(T-D) - g TanY (3.13) 

W V Cosy V 

h' = Tany (3.14) 


where 


I _ ^ 

dx 

Note that equation (1.17) corresponding to fuel consumption 
has been deleted in the above set. The Euler equation for 
the range problem is 




— 1- f/a - [T-D]) 

(T-D)U3h V av/ J 


(3.15) 


The initial value of the path angle y^ may be considered to 




be playing the role of X 
Jacobi-test matrix (3.12) becomes 


and since Xj^ is constant, the 


*av 

av " 


'av 

av ' 

ax. 

^0 

3^0 


^0 

ay 

o 




1 

0 


(3.16) 
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It is clear that the rank of the test matrix (3.16) is 

determined by the term . Therefore, if the sign of 

® o 

this quantity changes along the Euler solution, it is 
indicative of the existence of a conjugate point. 

Considering, next, the case of constant-density- 

atmosphere, the altitude becomes an ignorable state variable 
and the equations (3.13) through (3.15) becomes 


V’ = g(T-D) - g Tan? 
W V dosy V 

y = - 1 g a (T-D) 

(T-D) V W 


av 

To obtain analytical approximation to the term -=■ , the 

®o 

equations (3.17) and (3.18) are linearized 'about an 
equilibr^ium point defined by 


(3.17) 

(3.18) 


and, V obtained from 


resulting in 


Jf = Sin 


-1 


m 


a_(T-D) = 0 

av 


(3.19) 


(3.20) 


6V’ - - 6V - 


Sr = A, 5V 


(3.21) 

(3.22) 


where 


^1 " 


g IiD 

w Cosy V* 


“ 9 -Tany 


(3.23) 
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A, * - fg(T-D) SinY - g Sec^yl 
[ W V Cos*y T J 

(3.24) 


A, * - 9 _ 

a* (T-D) 


(3.25) 

V(T-D) 

av* 



upon further 

simplification. 

equations (3.23) 

through (3.25) 

become 

Al = 0 


(3.26) 


*2 = ^ 


(3.27) 


^2 ^ - yrr 

g (T-D) 

m T\\ <\ tt 2 

(3.28) 


Using Laplace transforms, equations (3.21) and (3.22) can be 
written as 


6V(s) = - A 2 6!T(s) 

*s” 

fiy(s) = ^ 6V(s) + 

s s 

Note that 6V^ =s 0. Combining (3.29) and (3.30), 


6V(s) = - 




Or 


5V(s) = - ■ A. 

s* + A-A- 
o 2 J 


(3.29) 

(3.30) 

(3.31) 


(3.32) 
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The roots of the denominator polynomial in expression (3.32) 
are either real or complex conjugates depending on the sign 
of the term A 2 A 2 • These two cases are considered separately 
in the following. 

( i ) ^ 2^3 ^ ® roots are + , real and symmetric 

about the imaginary axis in the complex plane. In this 
case^ inverse Laplace transform of (3.32) yields 

6 V - A, Sinh (Ox) 
o 

where 0 = >/-A^A^ • 

It is clear that the expression (3.33) is zero only at x = 
0. Consequently, if A 2 A 2 < 0, the rank of the test matrix 
(3.16) remains constant along the Euler solution and 
conjugate points will not occur. 


(ii) A 2 A 2 > 0, the roots are ± i . Inverse 

transformation of (3.32) then results in 


^ = “ A^ Sin (l-x) 

O 

where 'i' - ^ ^ 2 ^ 3 ' 


(3.34) 


In this case, the expression (3.34) changes sign at 


X 


n?T 

f 


(3.35) 
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and conjugate points will occur for sufficiently large range 

X. 

From the above, the condition for non-occurrence of 
conjugate points is 

A 2 A 3 < 0 <3.36) 

or 

_ (g/V)» 1 iL<T-D) < 0 (3.37) 

(T-D) 3^ 

For the glide problem, T - 0, and the expression (3.37) 
yields the following result : 

M 2 n 

if — > 0 , conjugate points will not occur. (3.38) 

The inequality (3.38) specifies a drag vs airspeed 
distribution which will be satisfied in all but very unusual 
aircraft configurations. This expression requires the 
aircraft to operate at the airspeed corresponding to minimum 
drag to maximize the range in gliding flight, and is 
consistent with engineering intuition. Note that the result 
(3.38) can be obtained by using a version of the energy 
model in conjunction with Mancill's theorem (Ref 90). 

Hence, the Euler solution for maximum-range glide under 
constant-density- atmosphere assumption satisfies the 
strengthened Jacobi necessary condition. Since these 
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trajectories also satisfy the strengthened Legendre-CIebsch 
necessary condition, the Euler solution affords a relative 
maximum for this problem. 

A similar exercise is next carried out for the time- 
range problem. 


3.2.2 Coni ugate-Po int Test for the Time-Range Problem 

The optimal-control problem in this case is the maximization 
of the final altitude for a specified range with fixed final 
time. The equations of motion (3.13) and (3.14) will be used 
in the following analysis also. With the interpretation of 
H as the time multiplier, the Jacobi test matrix for this 
problem becomes 


av av 

av 

ho o 


at at 

at 

ax, aH^ 

ho o 

'"o 

aXv. ax. 

ho o 

ax. 

Or 


av av 

av 

ho 0 


at at 

at 

ax. aH„ 

hn 0 

aTT 

o 


10 0 


(3.39) 


Hence, the sign of 
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av at - |v at (3.40) 

o o o o 

evaluated along the Euler solution determines the rank of 
the test matrix (3.39). If the expression (3.40) changes its 
sign along the Euler solution for the time-range problem, it 
is indicative of the existence of a conjugate point. 

Note that time appears in this problem as a state-like 
variable with 


t' = 1 (3.41) 
V Cosy 

As before, a prime on the variables denote differentiation 
with respect to the range variable x. The analysis given in 
the following will employ the constant-density-atmosphere 
assumption, permitting the deletion of expression (3.14) 
from the system. In this case, the Euler equation is 


y = Cosy 


V* (T-D) 


9 (H/X,) 3_[V(T-D)1 

^ av 


g ^(T-D) 
^(T-D) 3V 


(3.42) 


The equations (3.13), (3.41) and (3.42) are next linearized 
about an equilibrium point at a particular altitude, 


5V* = a^ 6V - a- 6y 
o 1 

6t' = - a2 6V + ag 6y 
6y* = a^ 6V - 


(3.43) 

(3.44) 


a- sy + ac 6H 

9 O O 


(3.45) 
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where 
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~ g ® (T-D) 

® W V Cosy av 


a, = 


= 


a- = 


3 . 

V 


V* Cosy 
Siny 
V cos*y 
- Cosy g 
V“ 

+ a_(T-D) 

av 


(H/X^) 


g 

V*(T-D> 


V(T-D)* 


V(T-D) av 


[‘ 

a_(T-D)(* 

av ( [ 

[ Co 


- Cosy (H/X ) 
V ^ 


- Cosy f H/X 
V 


+ ...g _ ^(T-D) 1 c^ (H/Xjj) 


1 

>'] 


= a. 


Sg = Cosy 


V^(T-D) X 


2 L-{T-D) + (T-D) 

X..1 av 


(3.46) 

(3.47) 

(3.48) 

(3.49) 


(3.50) 

(3.51) 

(3.52) 


Equations (3.43), (3.44) and (3.45) constitute a linear, 
constant-coefficient system which can be put in the 
following form using Laplace transforms (initial conditions 
on 6V and fit are zero) 


fiV(s) = 

«^o 


- a. 


s* + (a^a^ - a^ag) 


(3.53) 
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6V(s) = - a,ae 

fiHo — 

s[s“ + (aj^a^ - 

6t(s) = - [(a^ag - a^a^) - a^s] 

«o 

s[s* + (a^a^ - a^ag)] 

^(s) = - ((a a- - a^a,) - a-s] 

sMs* + (a^^a^ - 

putting = (a^^a^ - 

and 

T ss -ag 


(3.54) 


(3.55) 

(3.56) 


(3.57) 


(3.58) 


and cancelling out common constants in the niAmerator, the 
equations (3.53) through (3.56) can be brought to the form 

6V(s) = wj (3.59) 

XI 



, (3.60) 

6V(s) = 

««o - 

s(s* + ( 1 )^) 


5t(s) 

®"o 


(1 + Ts) u? 

n 


(3.61) 
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fit(s) = (1 + Ts) uj 



S*(S* + (0^) 


ORiGiNAL PAGE 5^ 

OF POOR QUALITY (3.62) 


Equations (3.61) and (3.62) may be further simplified using 
the expressions (3.59) and (3.60). 


6t(s) = 5V(s)+ T 6V(s) 

6ir^ 6H 

o o o 

6t(s) = u? " 


n 


s*(s* + w“) 


T 6V(s) 
5H 

o 


(3.63) 


(3.64) 


Equations (3.63) and (3.64) imply 


6t(x) 

«"o 

6t(x) 


6V(x) + T 6J7(x) 
6H_ fiy_ 


= L 


o 

-II 


u 


n 


(s* + w^)s*' 


+T ^(x) 
o 


using (3.65) and (3.66) in (3.40) 


- 6V 3t + 3V 3t = 


o o 
-6V(x) L 
«^o 


“li 


3H 3y 
o o 


n 


s^(s® + u^) 


+ l6y(x)f 


(3.65) 


(3.66) 


(3.67) 
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And consequently, one needs to obtain the inverse transform 
of only three transfer functions, namely 

§V{3) , §W{B) , 

6Z 6H " 

o o * 

s^(s* + w^) 


Two distinct cases can be seen to arise here, 


(i) > 0, the roots of the denominator polynomial is a 

complex conjugate pair and the inverse Laplace 
transformation yields 


av at - av at = 


ay 

o o 


ay 3H 
o o 


2 - X Sin(w^x) “ 2 Cos(w^x) 
n n n 


(3.68) 


The expression (3.68), after being zero at x = 0, will 


subsequently become zero at 


(3.69) 


Hence, in this case the rank of the Jacobi test matrix 
(3.39) will drop at the range defined by expression (3.69) 
implying the existence of a conjugate point 


(ii) < 0, the roots of the denominator polynomial are 
real, distinct and symmetric about the imaginary axis in the 
complex plane. In this case 


w jt - av at - 
3H ay ay an 

o o o o 
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2 + X Sinh(nx) - 2 Cosh(fix) 


where = 



(3.70) 


The expression (3.70) is zero only at x = 0. Consequently, 

the rank of the test matrix (3.39) will not change along the 

Euler solution and conjugate points will not occur. 

2 

The expression (3.57) for u is next evaluated with 

n 

the multiplier ratio H/X chosen from each of the three 

equilibrium points discussed in section 2.2. Unlike the 

maximum-range glide situation, the expression (3.57) is too 

involved to be amenable to analytical treatment. This 

2 

forces one to evaluate n numerically. It is found that 

n 

2 

is less than zero only in the third equilibrium regime, 
i.e., H/X„ values to the right of the V(T-D)__„ point. 
Hence, the oscillatory solutions shown in Figs 6 and 7 fail 
the conjugate point test for durations exceeding one period 
and furnish a minim\im only for short-duration problems. On 
the other hand, Euler solutions generated with the time- 
range multiplier ratio H/X chosen to the right of the 
airspeed for V(T-D) satisfy the strengthened Legendre- 

Iti aX 

Clebsch and Jacobi necessary conditions, and hence are 
optimal trajectories for the time-range problem. 
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3.3 CONCLUSIONS 

The Legendre-Clebsch and the Jacobi necessary conditions for 
the range optimal climb and glide, and the time-range 
optimal problems were examined in this chapter. It is found 
that the minimum-range climb or the "steepest climb" problem 
without any constraint on path angle possess no proper 
minimum or even a lower bound. If an artificial limit on 
path angle, say -90° < y < 90°, is imposed, the minimum- 
range trajectory consists of vertical up-down flight 
segments and any altitude-airspeed pair can be reached in 
zero range if there exists a positive margin of thrust over 
drag. On the other hand, the maximum-range glide path 
satisfies the strengthened Legendre-Clebsch and Jacobi 
necessary condition if a proper minimum-drag point exists. 

The Euler solutions for time-range problem meets the 
Legendre-Clebsch necessary condition with a margin. However, 
the strengthened Jacobi necessary condition is satisfied for 
"long" paths only by non-oscillatory solutions generated 
with the time-range multiplier ratio H/X chosen to the 

A 

right of airspeed corresponding to V(T-D) in Fig. 5. 

ni AA 

The analysis presented so far employed the constant- 
density-atmosphere assumption. In the next chapter, a 
numerical study of Euler equation using typical high 
performance aircraft (F-4) data is given. The discussions 
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there will also include a numerical conjugate-point test for 
the time-range optimal-control problem. 



Chapter IV 
NUMERICAL SOLUTIONS 


With the insaght gained for the range and time-range optimal 
problems with altitude dependence of thrust and drag 
supressed, a numerical study of the more general case in 
which the aerodynamic coefficients are functions of Mach 
number and the thrust is Mach- altitude dependent, is next 
undertaken. The data for a version of the F-4 aircraft with 
afterburner operative is used in this study. The thrust data 
is presented in Table. 1, and the zero-lift drag coefficient 
Cq^ and the induced drag coefficient are given as functions 
of Mach number in Tables 2 and 3 respectively. A cubic- 
splir.e representation (Ref 104) is used to interpolate the 
aerodynamic coefficients. The drag coefficient is computed 
as 



+ K 


Where Cj^ = 2 W/(p V^’ S) 

The drag is then obtained as the usual product of drag 
coefficient, dynamic pressure and the aircraft wing area. A 
cubic-spline lattice (Ref 104) is used to interpolate the 
value of thrust at a given altitude and Mach number. 
Atmospheric density and speed of sound as functions of 
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altitude are interpolated from standard atmosphere tables 
using cubic splines. The system differential eguanions 
(1.14) through (1.17) and the Euler equation (2.13) are 
integrated using a fifth-order Runge-Kutta-Verner method 
with variable step size. 

The level- flight envelope for the aircraft under 
consideration is shown in Fig. 10 along with the "energy 
climb" path/ generated by determining the maximum of 
specific excess power at each energy level. A discontinuity 
in the energy climb schedule due to transonic drag rise may 
be noted (Ref 44) . 

4.1 RANGE-OPTIMAL TRAJECTORIES 

Range-optimal trajectories are generated by setting the fuel 
multiplier and the variational Hamiltonian H to zero in 
the general Euler equation (2.13). If the initial data is 
given, the numerical solution for the system (1.14) through 
(1.17) and (2.13) can be generated. The full- thrust case is 
considered first. 

Typical time histories of altitude, airspeed and path- 
angle for the case in which a positive margin of thrust over 
drag exists, is shown in Figs 11, 12 and 13 respectively. As 
observed in Chapter II, the solutions are oscillatory, with 
the amplitude of oscillation depending on the particular 


65 


choice of initial conditions. For instance, if the initial 

value of altitude and airspeed are such that the partial 

derivative of (T-D) with respect to altitude, with specific 

energy held constant, vanishes, the oscillation has zero 

amplitude. As the initial conditions move away from the path 

^(T-D) 1 = 0, the amplitude of oscillation in 

3h 1 e = Constant 

the Euler solution increases. If one were to superimpose the 

numerical solution on the locus of (T-D) points in the 

max 

airspeed-altitude plane, the oscillations would be found to 
take place about this locus. It will be seen later in this 
chapter that the oscillatory paths do not satisfy the Jacobi 
test for durations exceeding half a period and hence, are 
not optimal. 

The situation with thrust set to zero is interesting, 

however. In this case, the Euler solutions are non- 

oscillatory and highly sensitive in nature. A maximum-range 

glide path is shown in the airspeed-altitude chart in Fig. 

14. This path is generated by adjusting the initial value of 

flight-path angle and integrating the system equations 

until the resulting trajectory passes through the desired 

terminal conditions. In the present case, the adjustment in 

y had to be carried out to 10-or 12 -decimal-Dlace accuracv 
o ‘ 

on an IBM-370/158 computer with quadruple precision. The 
remporal evolution of altitude, airspeed and flight-path 
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angle corresponding to this path are given in Figs. 15, 16 
and 17 respectively. These trajectories will be seen to 
satisfy the Jacobi necessary condition and are optimal for 
maximum- range glide. Note that jumps in path angle must be 
permitted at the initial and final points unless the 
specified boundary condition on happens to be the same as 
that emerging from computations. 


4.2 TIME -RANGE OPTIMAL PATHS 

The numerical investigation of time-range optimal paths is 
considerably more complex than for the range-optimal 
problem, mainly due to the involvment of the time-range 
multiplier ratio H/X , Of the various time-range optimal 

X 

paths, a trajectory of particular interest in practical 
applications is the climb-dash intercept path (Ref 13). 
This trajectory occurs in aircraft pursuit problems as a 
transient leading to the dash point or high speed point on 
the level-flight envelope, arising whenever the evader is at 
a sufficiently large distance from the interceptor. In this 
section, numerical computations carried out to generate the 
optimal climb-dash intercept trajectory will be discussed in 
detail. Note that the fuel multiplier X^ = 0. 

To obtain the value of the time-range multiplier ratio 
corresponding to the climb-dash path, an equilibrium 
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analysis along the lines of section 2.2 is first undertaken. 
Setting the time derivatives in equations (1.14) through 
(1,17) and (2.13) to zero, the time-range multiplier ratio 
H/X can be computed, with the specific energy E = h + 

A 

2 

V held constant. A sample plot resulting from these 
2g 

computations is shown in Fig. 18, with specific energy 
frozen at a typical value of 60,000 feet. Three regimes 
identified in section 2.2 can also be seen in this figure. 
The multiplier ratios, H/X , to the left of the point marked 
A (the this figure are all positive. 

Depending on the nature of thrust and drag, there may exist 
points at which H/X is zero. H/X values between point A 
and the point marked B ( the V(T-D)^^^ point) are negative. 
The time-range multiplier ratio H/X is positive to the 

A 

right of V(T-D) point. Numerical solutions with H/X 

in aX X 

values chosen to the left of V(T-D) point at various 

max 

energy levels indicated that they are oscillatory, similar 
in nature to the trajectories obtained in the full-thrust 
minimum-range climb problem. It will be seen subsequently 
that these oscillatory paths fail the Jacobi test for 
durations exceeding a period and furnish a minimum only for 
short-duration problems. 

The time-range optimal solutions generated with 
multiplier ratio H/X chosen to the right of V(T-D) point 
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at each energy level, are nonosci llatory and highly 
sensitive in nature. These paths satisfy the Jacobi 
necessary condition and are time-range optimal. To obtain 
the value of the multiplier ratio corresponding to the 
climb-dash path, a plot of the loc’..." of equilibrium points 
with a particular value of H/X , is made. The value of the 
multiplier ratio is adjusted next so that the locus of 
equilibrium points terminates at the dash point on the 
level-flight envelope. This value of H/X is used to 
generate optimal climb-dash trajectories. The equilibrium 
climb-dash path along with the "energy climb" schedule is 
shown in Fig. 19. The flight envelope is superimposed on 
this figure for clarity of presentation. As one decreases 

f 

the value of the range multiplier X it is clear that the 
equilibrium climb-dash path will approach the energy climb 
schedule. In the limit X =0, i.e., final value of range is 
no longer of interest, the two paths merge, a result 
consistent with engineering intuition. 

With the time-range multiplier chosen from the 
aforesaid analysis, all that remains to obtain the optimal 
trajectory is to determine the initial value of the control 
variable for a given set of initial conditions on 
altitude and airspeeo. An Euler solution for initial values 
of altitude and airspeed close to the equilibrium climb-dash 



69 


locus is shown in Fig. 20. The level flight envelope, energy 
climb schedule and the equilibrium climb-dash locus are also 
shown in this figure. It can be seen that the Euler solution 
stays close to the equilibrium locus. To determine this 
trajectory, an iteration was undertaken on the initial value 
of the control variable 2T. With quadruple precision on the 
IBM-370/156, had to be determined to 13 significant 

digits. To illustrate the sensitivity of the Euler solution 
to the initial value of the path angle Jf, the last digit of 
is perturbed in the positive and negative sense with 
trajectories 1 and 2 shown in Fig. 20 resulting. A few more 
Euler solutions with initial conditions far removed from the 
equilibrium locus are shown in Fig. 21. The temporal 
evolution of altitude, airspeed and path angle corresponding 
to the trajectory "a" in Fig. 21 is given in Figs. 22, 23 

and 24 respectively. 

An interesting feature of the Euler solutions for the 
Climb-dash intercept problem is that they funnel rapidly 
into a certain corridor in the airspeed-altitude chart, in 
the vicinity of the equilibrium locus corresponding to 
unaccelerated flight. This feature of the solution family 
may be exploited to simplify the computation of optimal 
trajectories . 
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4.3 NUMERICAL CONJUGATE -POINT TEST 

To carry out the conjugate-point test for the range-optimal 

and time- range optimal paths discussed in the earlier 

sections, the matrices (3.16) and (3.39) need to be 

evaluated along the Euler solutions. The numerical 

evaluation of the elements of these Jacobi test matrices 

evidently require computer codes considerably more complex 

than those required for the calculation of Euler soltions. 

As an alternative to the numerical solution of equations of 

variations, Cicala (Refs 92 and 93) suggested a scheme in 

which the partial derivatives with respect to the are 

■‘‘o 

calculated approximately in terms of difference quotients. 
Thus, small increments in the initial X^ are employed in the 
evaluation of neighboring solutions of the system (1.14) 
through (1.17) and (2.13). A computer code was assembled 
based on Cicala's suggestion (Ref 93) to evaluate the 
partial derivatives required in matrices (3.16) and (3.39). 
This code generates three trajectories corresponding to each 
multiplier, che first being the nominal and next two, the 
neighbor.i.ng trajectories obtained by perturbing the initial 
value of the multiplier in the positive and negative sense. 
The required partii'-I derivatives are then computed using a 
central differencing scheme (Ref 105). 
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One of the main difficulties encountered in this 
approximation is the errors arising from higher-order 
effects. Such errors can be controlled to a certain extent 
by checking the linearity of the x^(t) differences versus 
tile magnitude of the corresponding increment in the initial 
X value. This check has been incorporated in the computer 
program. 

The oscillatory minimum-range climb path is tested for 
conjugate points first. It is then found that the term 
3V/sy^ in the matrix (3,16) changes sign every half cycle of 
oscillation indicating the appearance of a conjugate-point. 
Hence these paths are non-minimizing. The maximum-range 
glide paths, however, satisfy the conjugate-point test. This 
result is consistent with that obtained in section 3.2.1 for 
the restricted case of constant-density-atmosphere. 

For the time-range problem, the conjugate point test is 
carried out for various values of H/X picked from the three 
equilibrium regimes at a particular energy level. It is then 
found that the oscillatory solutions obtained with H/X 
values from the first two equilibrium regimes indicate the 
existence of a conjugate-point after a cycle of oscillation. 
The non-oscillatory trajectories corresponding to H/X 
values on the right of V(T-D)^^^ at a particular energy 
level, on the other hand, satisfy the no-conjugate-point 
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condition. Hence, the climb-dash trajectories presented in 
Fig. 21 are optimal. 

4.4 CONCLUSION 

Numerical solution of the Euler equation for range-optimal 
and time-range-optimal problems was presented in this 
chapter. A method for choosing the time-range multiplier 
ratio H/X for the climb-dash mission was discussed in 
detail, A numerical conjugate-point test based on a scheme 
sxiggested by Cicala (Ref 93) was also presented. 

It is found that the maximum-range glide paths and non- 
oscillatory time-range solutions satisfy the Jacobi 
necessary condition and are optimal. Oscillatory time-range 
trajectories fail the Jacobi test for durations exceeding 
one cycle of oscillation. 

Numerical solutions of the Euler equation and numerical 
conjugate-point test has essentially reinforced the. 
conclusions arrived at in Chapters II and III, wherein the 
constant-density-atmosphere asssumption was used in the 
analysis , 



Chapter V 
CONCLUSIONS 

In this work, optimal flight in the vertical plane with a 
vehicle model intermediate in complexity between point-mass 
and energy models was studied. Flight-path angle takes on 
the role of control variable, and range-open problems 

feature subarcs of vertical flight and singular subarcs as 
previously reported. 

Minimum- range climb problem (the steepest climb of Ref 
5) has been found to have no minimum, not even a lower 
bound. In Ref 5, the steepest-climb problem was studied 

using the Green's theorem device of reference 4 and 90. 

There is an important difference in the vehicle modeling 

from that of the present work which should be noted as a key 
to resolving disparities between the character of optimal 
paths emerging. The analysis of Ref 4 and 5 in essense 

replace Cos3T in equation (1.16) with unity so that the 
problem solved is maximum altitude in a given distance (i.e. 
arc length) rather than in a given range. This is a 
necessity with the linear-integral approach which can 
accomodate only problems of dimension two and a very special 
form of state equations, A detailed discussion of the 
linear-integral approach and its relation to singulai 
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optimal control is given in an appendix. The solution to the 

distance-climb problem consists of a central path flown 

along a (T-D) locus in the airspeed-altitude plane with 
max 

vertical climb and dive tra.nsitions at the ends to meet 
specified boundary conditions. Imposition of a constraint 
on the path angle y -90^ -5 ^ 90°, showed that the 

"solution” to the ste.^pest climb can be constructed from 
vertical climb and dive flight segments alone, and that any 
altitude-airspeed pair can be reached in zero range. The 
possibility of this composite climb/dive path tending to a 
"chattering arc" was touched upon. Maximum-range glide 
paths emerging from the intex*mediate vehicle modeling were 
found to satisfy the strengthened Legendre-Clebsch and 
Jacobi necessary conditions. The use of constant-density- 
atmosphere assumption enabled detailed analysis of the Euler 
solutions and the Jacobi test for this problem. 

Since fixed- throttle minimum-fuel trajectories are not 
of interest in practical situations, main attention has been 
focussed on the time- range optimal control problem. For the 
special case in which thrust and drag depend only on 
airspeed, a plot of the ratio of time and range mutipliers 
H/X for equilibrium corresponding to unaccelerated flight, 
revealed the existence of three regimes. Positive values of 
H/X on the low-speed side of V{T-D) and all negative 
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values of H/X were shown to yield oscillatory solutions. 

X 

Although these meet the Legendre-Clebsch necessary 

conditions, they fail the conjugate-point test. Euler 

solutions with H/X chosen to the right of the V(T-D) 

X max 

point satisfy both the Legendre-Clebsch and Jacobi necessary 
conditions and are nonoscillatory in character. 

Numerical solution of the Euler equation and a 
numerical conjugate-point test for the F-4 aircraft data 
reinforced the conclusions arrived at in the analytical 
exercise . 

From a practical viewpoint, the time-range trajectories 
which terminate at the ”dash-point" on the level-flight 
envelope are of particular interest. The mutiplier ratio 
H/X corresponding to this point is determined using the 
locus of equilibrium points at each energy level 
corresponding to unaccelerated flight. With this value of 
H/X , Euler solution for any altitude-airspeed pair is 
obtained by iterating on the initial value of 2T, the control 
variable. Euler equations were obtained for various initial 
conditions. One observes that these rend to funnel rapidly 
into a certain corridor in the airspeed-altitude chart, in 
the vicinity of the equilibrium locus corresponding to 
unaccelerated flight. This feacure of the solution family 
can be exploited in practical situations to simplify the 
computation of optimal trajectories. 
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5.1 FUTURE WORK 

An important task to be undertaken is the comparison of the 
results obtained from the intermediate vehicle modeling with 
those generated by solving the complete two-point-boundary- 
value problem for aircraft. This will enable quantitative 
evaluation of the intermediate vehicle model and its 
usefulness in applications. 

Another interesting investigation would be the study of 
time-fuel-range problem with throttle control included. 
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Appendix A 


SINGULAR OPTIMAL CONTROL AND THE IDENTICALLY 
NON-REGULAR PROBLEM IN THE CALCULUS OF 
Vi^RIATIONS 


A. 1 INTRODUCTION 

In optimal-control problems featuring scalar control 
appearing linearly in the system differential equations, 
singular subarcs can sometimes arise. Along singular subarcs 
which are minimizing, the Generalized Legendre-Clebsch 
necessary condition should hold (Refs 72 and 88). A class of 
such optimal-control problems can be recast as identically 
non-regular problems in the classical Calculus of Variations 
if the dimension is low. Specifically, this transformation 
appears feasible if there are at most two non-ignorable 
state variables and one control variable. In general, the 
procedure involves a change in the independent variable 
under appropriate smoothness and monotonicity assumptions 
(Ref 94). (The phrase "classical Calculus of Variations" 
employed here refers to unconstrained problems, i.e., not to 
Lagrange-Mayer-Bolza problems.) 

For this class of problems, Mancill (Ref 90) has 
obtained conditions for a minimizing singular arc. In Ref 
90, Mancill made use of Green's theorem on line integrals to 
establish conditions for a strong relative minimum. Miele 
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{Refs 4 through 7 and 91) used the Green’s theorem approach 
for problems with control bounds, extended the technique to 
handle .'.soperimetric constraints and carried out 
applications to several flight problems. Haynes (Ref 95) 
discussed an extension of the Green's theorem approach to 
higher dimensions using exterior Calculus, for the optimal 
control of systems with n state variables and n-1 controls 
appearing linearly in the state equations. The question of 
existence of totally singular vector control is discussed in 
considerable detail in this work. Goh (Ref 96) examined the 
singular Bolza problem and noted the connection between 
Miele's work and the identically non-regular problem in the 
Calculus of variations. 

f 

This appendix deals with an evaluation of Mancill’s 
work and its relation to the Generalized Legendre-Clebsch 
necessary condition. A critique on the nature of 
transversality conditions for this class of problems is 
presented. Three illustrative examples are also given. 

A. 2 IDENTICALLY NON -REGULAR PROBLEM 

The identically non-regular problem with fixed endpoints in 
the Calculus of Variations (Refs 80 and 90) is the 
minimization of an inregral of the form 
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j = 


/ 


[P(t,x) + Q(t,x)x]dt 


with 


(A.l) 


x(t^) = x^ and x(t 2 ) = X 2 (A. 2) 
Note that 

[P{t,x) +Q(t,x)x]^^ = 0 (A. 3) 

XX 

It is known that the Euler's equation for this problem is 
either an identity or a finite equation (Refs 80,86,90 and 
97). If it is an identity, the integral is independent of 
the path joining two fixed points and no proper minimum 
exists. On the other hand, if it is a finite equation, the 
Euler's equation is satisfied only along certain paths which 
in general do not pass through the specified end points. 

These functionals are sometimes called "degenerate" 
because the Euler equation for such functionals is not a 
differential equation, but a finite equation without any 
derivatives of the unknown function (Ref 93). 

Two theorems by Mancill (Ref 90), yield the necessary 
and sufficient conditions for a strong local minimum in 
these problems. These are presented in the following. 
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THEOREM 1. If E ^2 is of class D' and minimizes the integral 
J in the class of admissible curves joining 1 and 2, where 
P(t,x) and Q(t,x) are of class in R, then 

g2n-lp/gx2n-l _ g2n>lg^g^g^2n-2 ^ 

3^^P/3x^^ i. 3^^Q/3t3x^^”^, (I) 

if 3^P/3x^ = 3^Q/3t3x^ k = 1,2/3, 2n-2, along arcs 

interior to R, including all isolated points in common with 
the boundary of R; 

3^P/3x“ ^ 3^Q/3t3x^"^ (Ig) 

if 3^P/3x^* = 3^Q/3t3x^ k = 1, 2 , 3 . . . . r-1 , along arcs Eg in 
common with the boundary of R. , 

Let (l') and (iL) represent conditions (I) and (I-,) 

O D 

respectively with the inequalities > replaced by the strict 
inequality >. This is a familiar notation in the classical 
Calculus of Variations and it will be employed in this work. 

The first part of (I) with n = 1, is the Euler's 

necessary condition for this problem. The inequality in (I) 
with n ■- 1 , is derived from the second variation. For n > 1 
the conditions (I) are obtained from higher variations. 

THEOREM 2, If P(t,x) and Q{t,x) are of class in R and 

the conditions (I') and (l!,) are satisfied along an 
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admissible curve joining 1 and 2 , then E ^2 furnishes a 
strong proper relative minimum for the integral J in the 
class of admissible curves joining 1 and 2. 


It is implied in Theorem 2 that the Euler equation is 
not an identity. This Theorem is proved using Green's 
theorem on line integrals. Mancill has given two additional 
theorems on the necessary and sufficient conditions for the 
identically non- regular problem with variable end points. 
However, the interpretation of these in the light of modern 
optimal-control theory points to the violation of the 
smoothness assumption essential to the results in Mancill' s 
work. A detailed discussion of this is presented in section 
A. 4. 

At this point, it is perhaps interesting to compare the 

results obtained by Mancill with those of Miele (Refs 4 

through 7 and 91). The first part of condition (I) in 

Theorem 1 with n = 1 is termed the "fundamental function" 

u(t,x) in Miele' s work. The inequality in (I) appears as a 

specification on the direction of traverse along the 

extremal. Similarly, the condition (1-,) of Mancill also 

o 

appears in Miele' s work as a specification on the direction 
of traverse along the boundary of the admissible region, 
applicable whenever the arcs interior to the admissible 
region are non-optimal. 
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A. 3 ^ PROBLEM IN ^ OPTIMAL-CONTROL FORMAT 

With a short development it will be shown that with n = 1, 
the inequality in (I) is the Generalized Legendre-Clebsch 
necessary condition for q = 1. 

Consider the optimal control problem 


Min 


/ 


lP(t,x) + Q(t,x)u]dt 


(A. 4) 


subject to the differential constraint x = u. 

It is apparent that this problem is equivalent to the 
identically non-regular problem in the Calculus of 
Variations. Note that the control u is unbounded. 

To proceed via the "modern" approach one defines the 
variational Hamiltonian 


H(X,x,t,u) = P(t,x) + Q(t,x)u + Xu (A. 5) 

and forms the adjoint equation 

^ “ (A-6) 

From the expression (A. 5) for H, one has that along a 
singular subarc 


= Q(t,x|t} ) + X(t) 


{A. 7) 
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Differentiating this with respect to time, substituting k - 
u and using (A. 6) for X,one finds 


d_(Hu] = Qt(t,x) - Pj^(t,x) (A. 8) 

dt 

Differentiating with respect to time again, while using x = 
u, leads to 

* <2tx - (A-9) 


Hence the Generalized Legendre-Clebsch necessary condition 
for first order singular arc is 


3 fe 

5Tr 


= 2tx - *'xx « ° 


that is 


(A. 10) 


P^^>-Stx (A -11) 

The inequality (A. 11) is the same as that in condition (I) 
of Mancill. 

One notes that the inequality (I) of Mancill for n > 1 
is not equivalent to the Generalized Legendre-Clebsch 
necessary condition with q > 1 but is something more 

general . 
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A. 4 TRANSFORMATION TO CANONICAL FORL 

To investigate the situations in which specified boundary 
conditions are off the path defined by the conditions (I), 
and the variable endpoint problem/ a transformation approach 
discussed in Ref 99 is next employed. The identically non- 
regular problem is first brought into the Mayer format : 

y = P(t/X) + Q(t/X)u (A. 12) 

X = u (A. 13) 

with t^, t^ , ~ ^2' ^^^2^ ~ ^2 ^ minimum of 

y(t 2 ) is sought with y(t^) = 0. 

Next a transformation of state variables will be 
performed so that the state system has a special form. The 
new state variables are z and x and the system is to have 
the control variable u appearing in only one of the state 
equations, the one for x. 

The system is 

z = P(t,x) + ^(t,x) (A. 14) 

at 

X = u {A. 15) 

and the choice of z leading to it is 

2 =y+R(t,x) (A. 16) 


where 
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X 

R(t,x) = - y*a(t,OdC (A. 17) 

(Refs 72 and 99). The end conditions are t^, x(t^) = x^^, 
x(t 2 ) = ^2 specified as before. The initial value of z is 
z(t^) = R(t^, x^) and a minimum of z(t 2 ) is sought. 

Since there are no bounds on the control u, it can 
behave impulsively and x(t) can jump. If the equation (A. 15) 
is discarded and a solution sought in the class of functions 
x(t) piecewise continuous, x becomes control-like (Refs 
72,99). At points t^ < t t 2 , X minimizes the right member 
of equation (A. 14). 

X = Arg min [P(t,x) + 3^ (t,x)] ' (A. 18) 

X 3t 

possibly exhibiting jump discontinuities in the interior of 
the interval depending on the nature of the time dependence 
of equation (A. 14). The variable x will generally jump at 
the initial and final times to satisfy the end conditions 
unless the value emerging from expression (A. 18) happens 
fortuitously to satisfy them. 

The situation with endpoint freedom is interesting. 
Consider for example, t^ and t 2 fixed as before, but x(t 2 > 
unspecified. To minimize y(t 2 ), x should jump at the final 
time t 2 to the value 

x(t, ) = Arg max R(t,, x) 

^ X 


(A. 19) 
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This seems to be the nearest thing to a transversality 
condition that one can have with x control-like. 

A. 5 ILLUSTRATIVE EXAMPLES 

To gain a better appreciation of Mancill's work/ three 
examples are given in the following. 

(1) Two elementary examples : 


(a) Min 


t 


f 


/ 


x^ dt 


t 

o 


, subject to X = u 


x(t^) = x^ and x(t^) = x^ specified. 

Since there are no bounds on the control/ the 
differential constraint is inactive. Hence/ the problem in 
classical Calculus of Variations format is 


Min 


/ 


x^ dt 


t 


o 


(A. 20) 


With the identification of 


P ( t, x) 2 x^ 


(A, 21) 
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Q(t,x) E 0 (A. 22) 

The necessary conditions of Mancill (Ref 90) become, 

2x = 0 (A. 23) 


and 


2 ^ 0 (A. 24) 

The sufficient condition 

2 > 0 (A. 25) 

is met in the strengthened form along the arc x = 0 and 
hence, the trajectory x = 0 affords a strong relative 

minimum. The result (A. 25) was obtained in Ref 72 via the 
Generalized Legendre-Clebsch necessary condition. 

If the initial and final conditions are off the x = 0 
path, jumps in x are required at the end points. Such 
motions have no effect on the performance index. 

The next example is chosen to illustrate the necessary 
conditions of Mancill for n > 1. 
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(b) Min 



subject to X = u 


The conditions x(t^) = x^ and x{t^) = x^ specified. Since 
there are no bounds on the control variable, the problem in 
the Calculus of Variations format is 



The necessary conditions for 
4x^ = 0 

Hence x = 0 is the extremal. 
12x^ = 0 
24x = 0 

24 0 


(A. 26) 

a minimum are 

(A. 27) 

Further, 

(A. 28) 
(A. 29) 
(A. 30) 


Note that the sufficient condition, (A. 30) with strengthened 
inequality, is met for n = 4. 

Just as in the previous example, jumps in x must be 
permitted at the endpoints if the specified conditions --re 
off the X = 0 path. 
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OF POOR QUALITY 


(2) Minimum-time aircraft climb 

Following Miele (Ref 4 through 7), a model of aircraft 
in symmetric flight under the assumptions of constant weight 
and small path-angle is 

V = g(( (T-D)/W| - Siny] (A. 31) 

h = V Siny (A. 32) 


Differential equations for range rate and fuel-flow rate 
have been dropped from the system, since they are ignorable 
in this problem. The optimal-control problem is the 
minimization of time required to fly from an initial (V,h) 
pair to a final (V,h) pair, viz. 

(V^,h^) 


Min 


/ 


dt 


(A. 33) 


(V. ,h. ) 
' 1 1 ' 


Changing the independent variable from time to altitude, 


V' = ^ = q(T-D) - g_ 
dh W V Siny V 


(V^,h^) 


/ 


dh 


V Siny 


(V,,h,) 


(A. 34) 


Min 


(A. 35) 
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Substituting next for Sin? in (A. 35) from (A. 34), the 
problem in the classical Calculus-of-Variations format is 


(V^,h.) 

Min f W 
J V(T- 

(V,.h,) 


+ W V dh 

D) g(T-D) 


(A. 36) 


In this development, the monotonicity of the altitude 
variable has been tacitly assumed. If desired, SinJf may be 
constrained by defining an admissible region in the V-h 
space as suggested in Ref 5; however, this fails outside the 
Manciil model. Employing conditions (I) in Theorem 1, the 
necessary conditions for a minimum for arcs interior to the 
admissible region, are 


a 

j w 1 

= i- 1 

j W 1 

(A. 37) 

3V 

|g{T-D)( 

3h 1 

1 V(T'D)( 



j w ( 

2^ if 

i _ w ( 

(A. 38) 

<1 

j g(T-D)j 

3h3V 

/ V(T-D)j 

The expressions 

(A. 37) 

and (A.3S) 

may be put in the 

following form 




3 

[V(T-D)] 

1 

= 0 

(A. 39) 

3h 

|E = Constant 

!l 

3h== 

[V(T-D)] , 

i 

< 0 

E = Constant 

(A. 40) 
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The sufficient condition for a strong relative minimum, 
then, ir;. 


3^ [V(T-D)] 
3h* 


(A. 41) 


E = Constant 

This result was obtained in Ref 72 using the Generalized 

Legendre-Clebsch necessary condition. The expression {A. 39) 

corresponds to stationary points of excess power V(T-D) 

' 2 

along contours of constant energy h + _V_. Inequality 

2g 

(A. 41) implies that the stationary points of excess power 
along constant energy contours must be maxima, a result in 
accord with engineering intuition. 

If the endpoints are off the path defined by (A. 39) 
jumps in airspeed and altitude must be permitted to meet the 
boundary condition. With bounds on control, on the other 
hand, operation at one of the control limits is indicated. 

The minimum- fuel climb problem may be handled in an 
analogous manner. 


A. 6 SMOOTHNESS DIFFICULTIES AND THEIR IMPACT 
In Mancill's paper (Ref 90), and in classical Calculus of 
Variations treatments generally, the function x(t), which 
appears along with its derivative, x(t), as an argument of 
the integrand, is assumed to possess a first derivative 
which is at least piecev/ise continuous. The various theorems 
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of Ref 90 do not apply to discontinuous solutions of the 
type examined in the preceding section. In the classical 
setting one would say that no minimum exists in the class of 
admissible functions, but only a lower bound. For this 
reason, the classical Calculus of Variations is sometimes 
dubbed the "naive" theory (Refs 86,100), Indeed the 
classical treatments (Refs 80,86 and 97) focus entirely on 
the degenerate case in which the integral is independent of 
the path. 

One is faced with the choice between extending the 
theory to admissible x(t) piecewise continuous, or the 
introduction of bounds on control u(t). Neither of these 
were done in Ref 90, and hence, this work is of limited 
applicability due to the explicit smoothness hypothesis. 
Another unwelcome complication in Mancill's work is the 
incorporarion of state-inequality constraints, a relic of 
his earlier work on this special type of problem (Ref 101), 
which do not alleviate the smoothness difficulties. 

Treatment of variational problems with x(t) piecewise- 
continuous only has been given by Young (Ref 102) and Krotov 
(Ref 98 and 103). Bounded-control problems approached by 
Green's theorem have been studied by Miele (Refs 4 through 7 
and 91 K In Ref 91, the control bounds are imposed by 
defining an admissible region in the state-space. 
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A. 7 CONCLUDING REMARKS 

Mancill's two Theorems given in the present work are of 
interest and seem to have been ahead of their time. For the 
narrow class of problems considered by Mancill, the 
inequality (I) with n = 1 is equivalent to the generalized 
Legendre-Clebsch condition. Perhaps equally important was 
Mancill's introduction of the Green's Theorem device fc ; the 
study of problems of small dimension. 
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Table. 2 


Zero-Lift Drag Coefficient Data for F-4 Aircraft 


WE I GH r 


MACH NO. 

^Do 

0.0 

0.0205 

0 . 2000 

0.0205 

0,6500 

0.0205 

0 . 7500 

0.0205 

0.8000 

0.0205 

0.8500 

0.0206 

0.9000 

0.0210 

0.9500 

0.02112 

1 . 0000 

0 . 032(1 

1 . 0500 

0.0359 

1 . 1000 

a .037‘1 

I . 2000 

0.03811 

1 . MOOO 

0.0365 

1 . 6000 

0.0386 

I . 8000 

0.0387 

2 . 0000 

0.0397 

2,2000 

0 , 0(103 

2 . 1«000 

0. 01103 

2.6000 

0.0403 


35UOO.O LBS WING AREA = 530 FT 
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Table. 3 . Induced Drag Coefficient Data for F-4 Aircraft 


MACH NO. 

K 

0.0 

0.1980 

0.2000 

0. 1960 

0.6500 

0. 1980 

0.7500 

0,2010 

0 . 8000 

0.2050 

0 . 8500 

0.2110 

0 . 9000 

0.2180 

0.9500 

0.2280 

^ . 0000 

0.2390 

1 . 0500 

0.2520 

1 . 1000 

0.2650 

1 . 2000 

0.2970 

t.ilOOO 

0.3680 

1 . 6000 

0.4530 

J.flOOO 

0.5440 

2 . 0000 

0.6400 

2 . 2000 

0.7200 

2 .' i 000 

0 . 7440 

2,6000 

0.7440 
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for the Range Problem 



Fig. 3 . Flight-Fath Angle vs. Airspeed in Powered Flight 


for the Range Problem 



30.00 


(T-D)/W 


50.00 , . .70.00 
AIRSPEED, Ft/S 


90.00 

» 10 ’ 


no. 00 130.00 


vs. Airspeed 


A Typical Parabolic 


Distribution 


OF poon QuAurf 



H/X^ V.. Airspeed for Equilibrium Flight (Parabolic 
(T-D)/w distribution) 

ax 


.a: 


.00 .„i90.00 110.00 

*10 


ical Solution for H/X in the 


gime 








‘30.00 50. 00 70 . 00 , 90 . 00 110.00 

AIRSPEED, Ft/S 3K I 0 


Fig. 8 . Representative Analytical Solution for H/X in the 

X 

Third Equilibrium Regime 
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Fig. 9 . A Parabolic distribution of Specific Excess Power 


vs. Airspeed 
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Fig. 10 , Level— Flight Envelope 'and the "Energy Climb" 

schedule for F-4 Aircraft 


A Energy Climb Path 
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M.OO 88.00 132,00. 17S.00 220.00 

VELOCITY.Ft/s «10‘ 

Flight Envelope, Energy Climb Schedule, 
Equilibrium Locus and a Climb-Dash Euler Solution 


A : Energy Climb Schedule 


B : Equilibrium Locus 


b : Climb-Dash Euler Solution 


2 : See P. 69 for details 







PIRSPEED.Ft/s xlO* 

150.00 164.00 178.00 192.00 206.00 220.00 






